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Abstract 

We propose a simple model of a waveguide network designed following 
the growth rule of a Vicsek fractal. We show, within the framework of real 
space renormalization group (RSRG) method, that such a design may lead 
to the appearance of unusual electromagnetic modes. Such modes exhibit 
an extended character in RSRG sense. However, they lead to a power law 
decay in the end-to-end transmission of light across such a network model 
as the size of the network increases. This, to our mind, may lead to an 
observation of power law localization of light in a fractal waveguide network. 
The general occurence of photonic band gaps and their change as a function 
of the parameters of the system are also discussed. 



PACS: 42.25.Bs, 42.82.Et, 61.44.Br, 72.15.Rn 

1 Introduction 

Propagation of electro-magnetic waves in dielectric materials has received extensive 
attention in recent years PP-[Zj. The existence of photonic band gaps (PBG) in the 
transmission spectra of such materials is of particular interest and is of importance 
in both fundamental science and technological applications. These studies constitute 
principal part of mesoscopic physics, a major focus being the possibility of localiza- 
tion of light jSlini, a direct evidence of which has already been reported for a strongly 
scattering media of semiconductor powders [Jj. The localization of classical waves 
is purely a result of multiple scattering in a random environment and free from the 
complications arising from interaction effects. Though, due to Rayleigh scattering, 
it is more difficult to localize classical waves compared to electrons. 

An alternative to the usual PBG systems may be obtained by forming network 
with slender waveguide tubes. This geometry does not need materials with high 
dielectric constants and yet, is capable of demonstrating Anderson localization of 
light. Zhang et al jH] have experimentally observed Anderson localization of light in 
a three dimensional network consisting of nearly one dimensional segments of waveg- 
uide. Vesseur et al P have examined the photonic band structure of a comb-like 
waveguide geometry. This work provides an example of a one dimensional photonic 
crystal enabling one to investigate the occurrence of localized states in such systems. 
The network system can produce large gaps even in one dimension, though it may 
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be sensitive to the structure of the unit cell. Recently, network models using serial 
loop structure have been studied in the context of propagation of electro-magnetic 
and acoustic waves [TUJUI]. Such theoretical investigations are additional justified 
as the present day nano-technology produce tailor-made waveguide networks with 
different geometries. 

Photonic band gaps also exist in the network systems without periodicity, such as 
fractal networks • Due to their self-similar structures and the absence of transla- 
tional invariance within the structures, fractals ^3 C3] are much more complicated 
to study compared to the periodic systems. Investigation of the propagation of clas- 
sical waves in a network model designed following a fractal geometry offers an in- 
teresting opportunity to observe the interplay of deterministic but non-translational 
order and interference effects. Unfortunately, this aspect has been grossly overlooked 
in the literature except for the recent work done by Li et al on a Sierpinski 
gasket network. 

In the present work, we propose a model waveguide network based on the geome- 
try of a Vicsek fractal ^3]- J7j. We study the transmission of electromagnetic waves 
through finite sized Vicsek Waveguide Network (VWN) using real space renormal- 
ization group (RSRG) and transfer matrix (TM) methods. A regular Vicsek fractal 
is already known to possess exotic spectral properties [THJ- [T7| . In addition to this, it 
has also been seen shown that an infinite Vicsek fractal can support a countable in- 
finity of extended wavefunction in the context of electronic transport [T%1 \W\ . Some 
of these states are 'anomalous' in the sense that the behaviour of the RSRG flow of 
the parameters of the system and the electronic transport are incompatible to each 
other JH]- This acts as an additional motivation for our work. We wish to exam- 
ine the appearence of photonic gaps in such a network system, the dependence of 
the gaps on the system size and the possible occurence of extended electromagnetic 
modes in such a deterministic fractal network. 

We find interesting spectal behaviour and formation of photonic gaps in a VWN. 
Additionally, we have been able to discern a number of eigenmodes in a VWN, some 
of which display a completely extended nature. For this latter set of states, the 
RSRG flow of the parameters gives a signature of extendedness, while the end-to- 
end light transmission decays following a power law as a function of the size of the 
network. This, to our mind, may lead to a possible power law trapping of light from 
an experimental point of view. We present out results below. 

2 The model and the method 

A 'unit cell' of a simple Vicsek geometry is obtained by considering a 'cross' with 
four dangling ends [T7j joined at a node, as a building unit. The next generation 
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Figure 1: Vicsek waveguide network in the (a) first and (b) second generation. The 
loops are attached to the open ends of a unit 'cross' as shown in (a). The patterned 
circles represent the nodes which are 'decimated' to obtain a Vicsek lattice of points. 

fractal is obtained by placing four such 'crosses' around a central cross and so on 
[T7| . A VWN may be constructed by considering each arm of a 'unit' cross to 
be consisting of a one dimensional waveguide segment. The dangling edges cause 
problems in realizing the boundary conditions from the standpoint of experiments. 
So, we attach loops of the same waveguide material at each end of a dangling arm. 
The first and the second generation VWN thus resemble what have been illustrated 
in Fig. H Each arm of a loop is of length a±. The other segments are of lengths a<i 
and 03 as shown in Fig. ^ 

Following Zhang and Sheng [20] the 'wave-function' ^,1+1 within any segment of 
length between the nodes i and i + 1 is given by, 



where, ipi is the amplitude at the zth node, and k is the absolute value of the wave 
vector in the loop. For an electromagnetic wave, k = where to and cq are the 
frequency and the speed (in vacuum) of the electromagnetic wave respectively, e 
is the relative permittivity of the dielectric medium, which we may assume to be 
real. A complex dielectric constant can easily be dealt with. The continuity of the 
wave function at the nodes and the flux conservation criterion are used to map the 
problem of wave propagation in such loops into an equivalent tight binding prob- 
lem of electron propagation on a one dimension lattice [20] -The resulting difference 
equation is, 




(1) 



(E - = t i>i+1 ll>i+l + ti-^i-! 



(2) 
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Here, we parametrize E = 2 cos klij+i and, define an effective on-site potential 



Si = 2cosM M+1 + £ 1 cot^ 1 ) >i + Y^cotO^l (3) 

m=l m=l 

and the effective nearest neighbour 'hopping integral' is given by, 

= I] . ( 4 ) 

m=i smt/ i i+1 

where, N^±i is the number of segments between the nodes i and z±l. 6^™^ = fcij'i+i 
is the corresponding phase acquired in the mth segment. The problem of wave 
propagation now becomes mathematically equivalent to the problem of transmission 
of an electron with energy E through the mapped lattice, where, — 2 < E < 2. 

It is now simple to eliminate the amplitudes ipi corresponding to the nodes de- 
noted by the patterned circle in [Fig. E(b)] in terms of the amplitudes at the re- 
maining sites. The result is a Vicsek lattice of points [Fig. 12(a)] . In the equivalent 
language of the electronic problem, we now have three distinct values for the on-site 
potentials, given by, 

4 / 1 

e A = 2 cos 6*3 + 2 cot 6»i + cot 2 + 



sin 2 6»i \E - 2 cos 3 - 2 cot 6 1 - cot 2 . 

ts = 2 cos 6*3 + cot 26*i + cot 9 2 

If 1 \ 

e c = 2 cos #3 + 4 cot 8 1 H ^— — (5) 

° d sin 2 fl 2 \E-2cos6 3 -2cot0i -cot0 2 / KJ 

and, the three nearest neighbour 'hopping integrals' are given by: 

2 / 1 



sin 9 1 sin 9 2 \E — 2 cos 03 — 2 cot 9\ — cot 02 
1 



sin0 



2 

t = -V ( 6 ) 

sm03 

where, T is the n.n. hopping integral between two neighbouring clusters. Here, 2a±, 
a 2 , a 3 are the lengths of the loop, the wire and the connector connecting the two 
neighbouring clusters respectively. Therefore, 9j = kaj with Oj = 01,02,03. This 
reconstruction facilates the use of RSRG to discern the extended and other modes of 
the lattice as well as to compute the transmittivity. The results are now discussed. 



3 Unusual eigenmodes of a VWN 

Let us consider a VWN at any arbitrary generation (converted into a Vicsek lattice, 
of points) clamped between two semi-infinite waveguide leads. Such a construction 
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e ^ o (a) ^ (b) (c) 



Figure 2: (a) A second generation VWN transformed into a Vicsek lattice of points 
on elimination of the patterned circles in Fig^ The effective (renormalized) on-site 
potentials and hopping integrals are shown, (b) The one-step renormalized lattice 
obtained from (a), (c) The final triatomic molecule to be placed between the leads. 



in its second generation is shown in Fig El Now, a VWN at any generation n + 2 
can be renormalized n times to obtain a second generation structure by decimating 
the appropriate vertices [T5] . The RSRG recursion relations for the on-site terms 
and the hopping integrals are given by, 

e A (n + l) = e A (n) + 



E - F 2 {n) 
t 2 2 (n) 



e B (n + l) = e B (n) , £ _ ^ 
e c (n + l) = e c (n)+4tl(n)F A (n) 
t x („ + l) - ^n)T(n)F l{ n) 



[E - e A (n))[E - F 2 {n)\ 
t 2 (n)t 2 l (n)T(n)F 1 (n) 



h{n + 1) " [E - e A (n)][E - F 2 (n)] 
T(n + 1) = T{n) (7) 



where, n denotes the stage of renormalization and 

E - e A (n) 



Fi(n) 



[E-e A {n)Y-T*{n) 



FM = e c (n)+ 2 * i(n) +t?(n)F 1 (n) 
h, — 6B\n) 



F 3 (n) 
FAn) 



E-e B (n) 



[E-e B (n)][E-ec(n)]-2%(n) 

E-e A (n)-tl(n)F 3 (n) 

[E - e A {n)\ [E - e A (n) - tl(n)F 3 (n)} - T 2 (n) 

It is seen that T remains fixed under RSRG iterations. 
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Figure 3: Distribution of amplitudes of the 'extended' wavefunction for the lead 
(dashed lines)-sample-lead systems corresponding to {a)E = e A — T and (b) E = 
e A + T. The transmission corresponding to these states across the sample decays 
with increasing system size as explained in the text. 

It is not difficult now to construct, by inspection, extended wave functions which 
have non trivial distribution of amplitudes along the X-axis. The wave vectors for 
these wavefunctions are obtained by setting 

E = e A (n)±T (8) 

As the RSRG process can proceed, in principle indefinitely, one can extract an 
infinity of different k- values by solving the equation (jSJ). A couple of such states are 
illustrated in Fig. El on a second stage fractal for E = e(n) ± T. It should be noted 
that €c in general flows to infinity with the progress of re-normalization. Therefore, 
it is necessary to choose the amplitudes at the central site of all the basic five-site 
cells to be zero in order to avoid the divergence of amplitudes at any arbitrary site. 
The 'extended' character of these eigenmodes are usually tested from a study of the 
flow of the 'hopping integrals' under successive RSRG steps. 

The recursion relations (J7J) reveal that if, E = e A (n) — T at any nth stage of 
re-normalization, then e A , cb, £i and t 2 immediately reach their fixed points i.e. 
e A (n + 1) = e A {n), e B (n + 1) = e B (n), tx(n + 1) = ti(n) and t 2 (n + 1) = t 2 (n) for 
all subsequent stages of renormalization. This implies that at any scale of length 
there is a nonzero overlap between the wave-functions at the nearest neighbour sites 
(at that scale), and that the amplitudes remain finite everywhere in even an infinite 
lattice. For example, solving E = e A — T we get ka = 0.547457038567240 and 
2.118253365362137 with a\ — a 2 — a 3 = a = 1. The hopping integrals t\ and t 2 
remain fixed at nonzero values (corresponding to their fixed points) in these two 
cases. Thus, from an RSRG point of view we expect an extended character of such 
states. To see whether such states are 'transparent' in regard of the incoming wave, 
we have to compute the transmission of an electromagnetic wave incident on a finite 
VFN through leads connected to the two-ends of the fractal lattice. 
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4 Transmission coefficient 



4.1 Formulation 

To calculate the transmission coefficient we first renormalize the Vicsek lattice n- 
times using the RSRG decimation method. The result is a Vicsek lattice consisting 
of five vertices and characterized by the n-times renormalized values of the 'on-site' 
potentials and the 'hopping integrals'. The recursion relations are given by Eq. (J7Jl. 
This renormalized version [Fig. EJb)] is now reduced to a 3-site cluster [Fig. |2Jc)] 
and the effective on-site potential of the central site becomes 

This triatomic 'molecule' is now clamped between two semi-infinite perfect waveg- 
uides, which are represented, following the discretization process outlined in Ref 
by a sequence of identical site potentials 



e = 2 cos 6 3 + 2 cot 6 3 (10) 
and nearest neighbour hopping integral 

to = -V ( n ) 
sin 3 

The decimation procedure is illustrated in Fig. |21 for a second generation fractal. 
The equations are easily obtained by artificially considering the lead to be consisting 
of identical materials of length 03 joined end-to-end. The transmission coefficient r 
is now given by |21j . 



[Pi2-P2i + (Pii-P22)cos6tf + (P 11 + P 22 ysm 2 6 3 1 ' 



where, 



P 
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P 
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[E - e A (n)} 2 [E - ~e c (n)} 2[E - e A (n)] 
tj{n)t t 
[E-e A {n)][E-~t c {n)} 



ti(n) 
P21 = — P12 

t [E-e c {n)} 



P 



22 



t\{n) 



4.2 Results 



(a) Let us analyze the wavefunction which corresponds to an extendedd state of 
the lead-VWN-lead system. For the values of wave vector correspond to the energy 
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Figure 4: (a) r vs ka 3 with a\ = 0.01, a 2 = a 3 = 1. (b) a x = 0.05, a 2 = a 3 = 1. The 
solid and the dashed lines in each figure indicate results for the second and the first 
generation VWN respectively. 
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E = €a — T, for example, we can use the recursion relations (J7J) to write explicitly 
the renormalized values of ec in the successive generations as, 



e c (l) = 5e c - cr 
e c (2) = 5 2 e c — ha — a 

and so on. Here 

is constant because €a, £b, fa, fa and T attain their fixed point values for n > 1. 
Continuing in this way we get, after n steps of renormalization, 

e c (n) = 5"e c - P„_i(a) 

where, 

P n -i(a) = 5"-V + 5 n -V + + a 

and 5 n = AT, the system size of the lattice at its nth generation. Therefore, we see 
that for large values of n, i.e., for large size of the original lattice ec(n) ~ iV as €a, 
6b, fa and t-i are invariant at each step of renormalization. So for large n, its simple 
to work out from Eq. (jHJ) that r ~ i indicating power law decay as a function 
of the total number of matrices. This is incompatible with the invariance of the 
hopping integrals under RSRG. The growth of ec(n) under RSRG implies that an 
electromagnetic wave travelling through the perfect waveguide leads and entering 
the fractal structure will experience larger and larger effective 'potentials' at the 
centres of the five site clusters. For large enough VWN systems we thus have the 
possibility of a power law localization of light. 



The existence of the dangling side units consisting of a wire and a loop which 
themselves form a fractal geometry, causes destructive interference between the 
waves propagating and reflecting back and forth along these side units. As for 
all these wave vectors the amplitude of the wave function is pinned at the value 
zero at the centerers of five site clusters at different length scale, standing waves 
are formed along the branches with nodes at the centers of one or the other five 
site cluster, which should localize the light wave and reduce the transmission. We 
call such modes 'unusual', and emphasize that these are eigenmodes not for a finite 
cluster, but for the lead- VWN- lead system, irrespective of its size. 

(b) We have also examined the features of the transmission spectrum of finite 
VWN's of arbitrary size. The appearance of transmission maxima separated by pho- 
tonic gaps is a general phenomenon. The fragmentation in the spectrum increasing 
with increasing size of the system. Figures (@J) and (0) exhibit the main character- 
istics where, in addition to the variation of r against the wave vector, we study the 
effect of the change of the size of the loops (cti). Three cases have been explicitly 
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Figure 5: r vs ka^ with a\ = 2, a 2 = «3 = 1. The solid and the dashed lines now 
represent the results for the first and the second generation VWN respectively. The 
appearance of new transmission zones is clearly seen between kas = and 2, 4 and 
6. 

shown with the size of the loop (ai) increasing from a very small value to a value 
comparable to the lengths of the other waveguide segments (a 2 and a 3 ). 

Though in every case, the exact distribution of values of the tau is sensitive to 
the individual values of the parameters, still some interesting behaviour is observed 
in certain cases. For example, setting a 2 = = 1 and a\ = 0.01 in arbitrary units, 
the transmission spectrum shows the presence of clusters of high transmittivity 
flanked by sharp delta-like peaks, even for a first generation VWN [dashed line in 
Fig. Ufa)]. The same set of parameters leads to a fragmented band structure as we 
move over to the second generation [vertical line in Fig. Efb)]. The spectrum reveals 
almost similar features as a\ is gradually increased from 0.01 to higher values. For 
a\ = 0.5, we get a clear signature of wider photonic gaps separating four clusters of 
high transmittivity. This is seen in Fig. EJb), where the dashed and the solid lines 
again correspond to the first and second stage VWN. 

When the dimension of the loop starts getting larger than the connecting waveg- 
uide segments, we find that new clusters of transmission maxima start appearing in 
the major photonic gaps of the earlier generation [Fig. The three major pho- 
tonic gaps corresponding to the first generation fractal get populated by two major 
photonic bands (dashed line in Fig. |3j) when the length of the loop exceeds the other 
two segments a 2 and a 3 . 

To make the appearance of the photonic gaps more pronounced, in Fig. |U] we 
show the variation of transmission coefficient r of a VWN in the first (dashed curve) 
and the second (solid vertical lines) respectively. The more or less fixed zones of 
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Figure 6: r vs a\ with a 2 = a 3 = 1 and k = n/2. 

zero transmission around a± — 1 and ai = 3 mark the broad photonic band gaps in 
a VWN. 

5 Conclusion 

We have shown that light transmission in a Vicsek fractal geometry is marked by the 
appearance of unusual extended-like eigenmodes for which a non-trivial distribution 
of finite amplitudes can be obtained in a Vicsek waveguide network clamped between 
semi-infinite leads. The transmittivity of such a lead-network -lead system shows a 
power law decay as the network grows bigger. Broad photonic gaps are obtained by 
appropriately choosing the network parameters. 
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